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The electronic states of a parabolic quantum dot in a magnetic field are studied with the inclusion 
of the spin-orbit interaction. The analitycal formulae for the ground state energy of the interacting 
\Q ' system are derived. The spin-orbit interaction is shown to introduce new features to the far infra- 

red absorption spectrum, where it leads to the splitting of the two principal modes. The results 
are compared with the charging experiments by Ashoori et al. and the far infra-red absorption 
measurements by Demel et al. 
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I. INTRODUCTION 



Recent progress in the semiconductor nanostructure technology has allowed for creation of the quasi-zero- 
■ dimensional (OD) systems known as quanturn-dptsa. These structures can be obtained e.g. by applying spatially 
[ modulated electric field to the 2D electron gasacl, or by embedding a small piece-,oi one semiconductor into another, 
\Q . with hieiier conduction band energy (this can be achieved e.g. through ctchinguLl, intcrdiffusionlj cp-sclf-organized 
growthBt2l) . The resulting lateral potential is often very well approximated by the isotropic parabolaliIliL2l. The con- 
finement in all three spatial directions leads to a discrete energy spectrum of the system jy&h characteristic excitation 
energies of the order of few meV, which can be studAcd using the spectroscopy methodsBtJ. 

Due to the so called generalized Kohn theoremolij for the parabolic confinement the electron-electron Coulomb 
interaction does not affect the resonance energy spectrum under the far infra-red (FIR) radiation. The FIR resonant 
spectrum of the correlated many-electron dot is therefore identical to that of a single electron and consists of two 
. degenerate modes. The degeneracy can be removed by a perpendicular magnetic field under which the two modes 
evolve into the inter-Landau-level (cyclotron) and intra-Landau-level frequencies. However, a number |Ofj-experiments 
seem to show a slight dependence of the two resonance-energies on the number of electrons in the dotou. Moreover, 
O ■ a small splitting of the two principal modes is observedo what reveals their additional sub-level structure. 

The theoretical investigations of the ground state of quantum dots containing 2-15 electrons were reported e.g. in 
Refs. 15 -2(]. In this paper we calculate analitically the ground state energy of a system of a larger number of confined 



electrons (~ 15-100) within a Hartree-Fock approximation. The spin-orbit interaction is included here in the manner 
analogous to that used for many-electron atoms and not via_the bulk band-structure parameters, since the considered 
system is strongly localized (diameter ~20 lattice Constantsa) and the small piece of crystal, to which the motion of 
electrons is limited, cannot be treated as an infinite periodical lattice. Instead, we propose to extract the spin-orbit 
coupling from the band-structure description and include it later on the level of the Bloch envelope wavefunctions 
by adding the appropriate term to the hamiltonian. As demonstrated by Darnhofer and R6sslcrE2l for 2 electrons, 
inclusion of the spin-orbit interaction through the band-structure parameters in the InSb dot leads to the similar 
effects on the electronic structure as obtained here for the GaAs dots (spin-orbit coupling in bulk InSb significantly 
exceeds that in bulk GaAs). 

The results obtained within proposed here framework seem to explain a number of experimentally .observed effects, 
like an appearance of the higher modes in FIR absorption and their anti-crossing in a magnetic fieldo, or the charac- 
teristic bumps in the magnetic-field dependence of the ground state energyQ. Let us also underline that we managed 
to fit very well the characteristic magnetic field at which the anti-crossing in FIR spectrum and the bumps in energy 
occur. The fact that this w.as impossible without including the spin-orbit interaction (see the numerical results and 
discussion by Palacios et aiEl for up to 15 electrons) seems to prove the importance of this effect. 
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II. MODEL 



The interaction between the spin and orbital angular momentum: er and 1, of an electron confined in a quasi 
two-dimensional quantum dot is included in the way analogous to that used in many-electron atoms, i.e. via the 
single-particle potential: 

V LS = a\ ■ a. (1) 

The coupling constant a is connected with the average self-consistent field < SU > acting on the electron, via the 
relation: 

a = (3 < SU > . (2) 
For a Z-electron atom the dimensionless parameter j3 is: 

Ze 2 \ 2 _ ( Z_ 



In the case of a quantum dot (3 will be treated as a fitting parameter, while the magnitude of the field < SU > will 
be estimated according to the electronic structure of the dot. 

Let us write the total hamiltonian of a system of many electrons, in the effective mass approximation, including the 
kinetic energy in the perpendicular magnetic field B, parabolic confinement of the characteristic frequency u>o, spin- 
orbit coupling as above, the Zeeman splitting for the effective ^-factor and the electron-electron Coulomb interaction 
controlled by the dielectric constant e: 



e 2 



i i^3 J 

In the above m* is the effective mass, r are the positions, p = —ifiV are the momenta, A = \B{y, — x, 0) are the 
vector potentials in the symmetric gauge, I and a are the projections of the orbital angular momentum 1 and spin er 
across the plane of motion. 

In the Hartree-Fock (HF) approximation the equation for the HF wavefunctions ip reads: 



[H B + Vi]i/>i(Ta) +J2 j dv ' A;(rcr,rV)^(rV) = e^ra), 



(5) 



where Hb is the hamiltonian of a single (non- interacting) electron for the field B defined in Eq.(Q), Vi denotes the 
Hartree potential: 

Vi = e -[ d r>^l, (6) 
e J |r'-r| 

with: 
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and A, is the Fock correction: 



A i (w,ro / ) = --U^ | r/ _ r i ■ (8) 

3 

Introducing the exchange operator G as: 

G^(ra) = E/ dr ' ^(rcsrVOV^iV), (9) 
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the HF equations (Eq.^J) can be written in a compact form: 

[H B +Vi + Gi] ^(ra) = e^ra). (10) 

Following the work by Shikin et aI0 we employ here the approximate formula for the charge density in the parabolic 
dot (of confining frequency loq), obtained within the classical regime and therefore applicable for large numbers of 
electrons N: 



l(r) f„(0)i^Ff„,., : 



< R 
for r > R 

with the charge density in the center n(0) = 3N/2irR 2 . The dot radius R, for the classical system given explicitely in 
Rcf. ^l] will be calculated here in a variational manner (from the minimum energy condition) in order to account for 
the quantum corrections. Using the formula ( |TT| ) the Hartree potential is calculated: 

3nNe 2 ( r 2 \ 
It It \ 

(we take this form also for r > R). 



V R (r) = 1 r (12) 
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III. QUANTUM DOT AT ZERO MAGNETIC FIELD 
A. Single-electron states 

Let us first consider the case of zero magnetic field. Neglecting for the moment the exchange term (which will be 
included later as a perturbation) we obtain the following Hartree equations: 

A + ~mXV + cda) ifc(r<r) = (e« - M™), (13) 



2m* 2 u j r ' v ' V 4eR 

where we use the effective frequency, renormalized by the Hartree term (Eq.[l2"|): 

% 2 =c 2 - *" Ne \. (14) 
4em*i? 3 V ; 

Eq.(|l3|) can be solved analitically and the obtained eigen-states are: 

ipi(ra) = ip nma (Ta) = <j) m {0)R nm {r)xa, (15) 

with the spin eigen-function \ a (with eigen- values a — ±h/2), the angular wavefunction of the angular momentum 
eigen- value m: 

cP m (e) = ^=e lmQ (16) 



/2tt 

and the orbital wavefunction: 



*-(r) = ^ jT^bm ( W ) ' e-^LtJ ( ' W ) . (17) 



V2 
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In the above L ™ are the Laguerre polynomials: 



1 H nr 

L^(z) = *"I"V (^+1^-*), (18) 
m! dz™ 1 - 



Zq = y/h/mQ,' is the characteristic length, n = 0,1,2,... is the principal quantum number, m is the azimuthal 
quantum number (|m| < n and the parity of m is the same as that of n), and n r — "~^ T "I is the radial quantum 
number. 

The eigen-energies associated with the eigen-functions ipnmo- are: 



3 



hQ (n + 1) + ama 



SirNe 2 
AeR ' 



(19) 



In the absence of the spin-orbit interaction (a = 0) they form degenerate shells labeled by n. Non-vanishing a splits 
these shells into the doubly degenerate sub-levels. 

Often used is a complementary (Fock-Darwin) representation: 



(20) 



with n± = 0, 1, 2, . . . The two sets of quantum numbers: [n, m] and rt_] are connected by the simple relations: 
n = ra+ + n_ and m = n+ — n~ . The orbital part of tpi is defined as: 



* (T)= 1 (° + ) m+ ( 6+ )"- c -'/2,g 



(21) 



where the raising operators a + and b + are: 



x + iy 



„ . d .d 

\dx +l dy 



x-iy , 



d d 



In hr- -i^r- 



l'o ° \dx dy 



(22) 



The eigen-energies labeled by n± and a are: 



1 N , 1, 3irNe 2 
£+ (n + + -)+ e _(n_ + -) + ^ r , 



(23) 



where e± = hQ' ± acr. 



B. Many-electron ground state 

The ground state energy of the system in terms of the N lowest Hartree eigen-energies £j, with i standing for the 
composite index [n, m, a], reads: 

where the subtracted integral represents the direct Coulomb energy of the system, counted twice in the summation 
of the Hartree energies £j. Introducing the Fermi energy sp separating the occupied and unoccupied Hartree energy 
levels in the ground state, and calculating the self-interaction integral we arrive at the formula: 

Eo= Y / Q(e F -s i )s i -^f. (25) 

i 

In the above is the Heaviside function. The Fermi energy is determined by imposing the fixed number of electrons 
N: 

N = J2 ®( £ f -£<)• (26) 

i 

The details of calculating the Hartree energy £q are given in Appendix and here we shall only present the final 
result: 



_ 9^V 2 e 2 2 3/2 / f3*N 

The radius of the dot i? can now be determined from the minimum condition: d£o/dR — 0, equivalent to the 
equation: 
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(28) 



where ag — th 2 /m*e 2 is the effective Bohr radius. Confining our considerations to the case of the large number of 
electrons, we can solve the above equation perturbatively with respect to the small parameter clb/R <§C 1. The zeroth 
order approximation Rq can be written in the form: 

» 3nNe 2 

R a = 1 9 29 ) 

and coincides with the classical result obtained in Ref. [2l|. Assuming the first order approximation in the form: 
R = i?o(l + 6), the correction 6 reads: 

-^('-^)- 

Disregarding higher corrections (non- linear in S), the effective confining frequency Q' can be found based on its 
definition (U4f): 

n' 2 = n 2 (i-^f], (3i) 



3G 



where we use the notation: 



2 _ lOOoB 2 

na ~27^R °- { } 

Finally, we arrive at the formula for the ground state energy in the Hartree approximation: 

£o = ^ + >o(l-^V /2 - (33) 



20ei? 3 u V 36 

The first term in the above equation, calculated in Ref. [2l|, is the classical energy of N interacting electrons confined 
in a parabolic well: 

9irN 2 e 2 t ,. 1 , 3J e 2 t f , , n{r)n(r') 

= / drn(r) ■ -m*ujf,r 2 H / dr / dr' , \' . (34) 

20ei? J V ' 2 2eJ J |r - r'| V ' 

The second term is the quantum correction and splits into the energy of oscillation with the frequency f2o : 

^hn Q N^ 2 = ©(£F - Ei) \m*nl < t\r 2 \t > (35) 

i 

and the spin-orbit interaction term. 

Let us now include perturbatively the exchange interaction, neglected so far in the Hartree approximation. As the 
first order correction A£ we shall calculate the average value of the exchange operator G, defined by Eq.(^), in the 
Hartree ground state obtained without including the spin-orbit interaction: 

A£ = - Si) < i\G l \i >| /3=0 . (36) 

i 

We have verified that the effect due to the spin-orbit interaction is indeed negligible here: for (3 — 0.3 the correction 
to A£ does not exceed 0.15 meV per electron for N < 100 (compare with the energy scale in Figjl]). Similarly, the 
second order correction proved to be smaller than the first order correction by the orders of magnitude (reaching 
merely 10 -2 meV), which is then a good approximation to the actual exchange energy. 
As shown in Appendix |^, Eq.(|3^) can be conveniently written as: 

/ 3 N w 

where <5 = 8{(3 = 0). 

Thus we have obtained the total ground state energy £ of the system of iV electrons confined in a parabolic well, 
including the kinetic energy, the direct and exchange Coulomb interaction, and the spin-orbit coupling: 

£ = £n + AS. (38) 

In Fig.|l| we present the average ground state energy per electron e = S/N plotted as a function of the number of 
electrons N. The two curves corresponding to the parjasneter /3 equal 0.3 and 0.6 are shown to bera reasonable 
interpolation between the classical result by Shikin ct alEil and the experimental data by Ashoori et alB. 



f 3 \ A 7 /V , 
A£ = -j= 1 1=) 5— (l-<5o), 37 
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C. Far infra-red absorption 



Let us now consider the selection rules for the optical transitions of the system under the far infra-red (FIR) 
radiation. Absorption of the FIR light, leading to the excitation of the electron droplet, has been a powerful tool in 
the experimental studies of quantum dotscTEl. 

Since the wavelength of the FIR light is much larger that the radius of the dot, one can use the dipole approximation 
for describing the interaction between the light and electrons. The probability of the optical transition between the 
initial (i) and final (/) states is proportional to the squared matrix element of this interaction: 

d% ~ | < /|eE • J2 - e<) r,|z > | 2 , (39) 

i 

where E is the electric field, uniform over the volume of the dot. The dipole matrix element dfi vanishes unless there 
is a pair of the HF states, one in the initial and the other in the final many-electron state, with equal spins and 
different by unity in each of the orbital quantum numbers [n,m]: 

cr f = <j\ \n f -n i \ = l, \m f -m i \ = l, (40) 

with all other corresponding HF states equal in the initial and final state. In other words, the absorption of a FIR 
photon leads to the excitation of a single electron from its (initial) HF state to another (final) HF state with the 
same spin a and the orbital quantum numbers changed according to Eq. ( |io| ) . Translating these selection rules to the 
Fock-Darwin representation we have: 

a-f = a 1 , n+ = n\ ± 1, n_ = n l _, 
or : a* = a' 1 , n S _=n l _± 1, n+ = n\, (41) 

i.e. the excited electron changes one of its orbital quantum numbers [n + ,n_] by unity. 
The above selection rules lead to the splitting of the resonance energy: 

£f - £ l = e ± = %n' ± X -ol, (42) 



The magnitude of this splitting a depends on the number of electrons according to Eq.(B5) 



IV. QUANTUM DOT IN A MAGNETIC FIELD 
A. Single-electron states 

Sketched in the previous section for the case of zero magnetic field the procedure of minimization of the Hartree 
energy with respect to the dot radius R, with later perturbative inclusion of the exchange interaction, has been also 
carried out for nonzero fields. The explicit form of the Hartree equations including the presence of a perpendicular 
magnetic field is: 

A^A + X -m* (o! 2 + ^ r 2 - huj - gf i B <jB + ala^j ^(m) - ^ - |^|y) ( 43 ) 
where lo c = eB/m*c is the effective cyclotron frequency and the zero field radius R appearing in the definition of 



Q,' B (14) is now replaced by R(B). We shall also denote the total confining frequency by: fi' 2 = fi'J + cj 2 /4, and its 



corresponding characteristic length by: I = yjh/m*£l' 1 . 

The eigen- functions of Eq.(|43]) are of the same form as given by Eqs.([l5]) and (|2^), only with the characteristic 
length replaced by I. The corresponding eigen-energies read: 

1 3nNe 2 
£nma = W. (u + 1) - -huj c m - g(i B <?B + ama + |^7gj > (44) 

or, in the other representation: 

/ 1 ^ , 1 s n 37riVe 2 

£ n+ n-a = £+{n+ + -) +£_(n_ + -) - gii B oB + ^ R , (45) 



G 



with e± = Q' ±(hu) c /2 + a<j). Since the Zeeman splitting is rather small for GaAs (g ~ \ yielding g[iB ~ 0.05 meV/T), 
which is the most common material used for the quantum dots, we shall neglect it in the further considerations. 

Including the magnetic field leads to the possibility of crossings between different energy levels £,*. Whether the two 
close levels e\ and £2 actually cross, or their crossing is forbidden, depends on vanishing of the off-diagonal matrix 
element of the operator describing the change of the hamiltonian due to small change of the field. Thus the condition 
for the allowed level-crossing is: 

r) J-f 

<1| — 12>^0. (46) 

The operator dH/dB commutes with the spin and inversion (r — > — r) operators. Its commutation with the angular 
momentum requires assumed here circular symmetry of the confining potential. Hence, while for the states of non- 
equal quantum numbers n and a we have the condition (|46| ) guaranteed, for a pair of states different only in m in 
general it is no longer true. This leads to the anti-crossing of levels, which can be taken into account by changing the 
formulae for eigen-energies: 

£nma = W! (ft + I) + m 

Analogously we have to modify the definition of a pair of energies e+ and £_, appearing in Eq.(f45|), into: e± = 
Q' ± \hu) c /2 + aa\. Let us underline that this rearrangement of levels is a perturbative approximation, beyond the 
Hartree-Fock approach. 



+ 
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(47) 



B. Many-electron ground state 



The many-electron Hartree ground state in the magnetic field £q(B) is defined analogously as in Eq.(^ 
the procedure sketched in Appendix ^ one can bring it to the form: 



£o(B) = 



20ei? n 



■-N 3 ' 2 h 



uj 2 u b 



3nNe 2 



1/2 



I - 



36 : 



Using 



(48) 



with the functions fs and ub defined in Appendix || by Eq.( DI6 ) and Eq.( DI 7 ), respectively. 

Analogously as for the zero magnetic field, the ground state radius of the dot R(B) can be found from the minimum 
energy condition: 



d£ (B) 
dR(B) 







which resolves into the equation: 



3iTNe 2 



ui u B = 



Aem*R(B) 3 
In the zeroth order approximation we obtain: 



I00a B 
27irR{B) 



1 - 



36 



Ro(B) 3 



3TrNe 2 



1 R3 



where Rq is given by Eq.(E 



iem*uj 2 UB 

The first order correction 5(B) defined as: R(B) = Rq(B)(1 + 6(B)) reads: 



6(B) = 



IOOqb 
81nR (B) 



1 - 



P 2 f 2 B N 
36 



(49) 



(50) 



(51) 



(52) 



In Figj^ we have drawn the radius as a function of the field. The dependence is fairly weak. When the magnetic field 
is increased, at low fields the radius also increases, and later, at higher fields, slightly falls down. The initial increase 
of the radius in the rising field is an effect due to the spin-orbit interaction, and vanishes for /3 = 0. 
Finally, the ground state energy in the Hartree approximation can be found in the form: 



£o(B) 



9nN 2 e 



2^2 



f 



2QeR (B) 3 



u 2 g 3 hQo ( 1 



36 



(3 2 f 2 N 



N 3/2 . 



(53) 
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Calculating the correction due to the exchange energy in the obtained above Hartree ground state is far more 
complicated for non-zero magnetic fields. Therefore the hypothesis is used, according to whici-|the kind of dependence 
of the exchange energy on the number of particles is not affected by the presence of the fieldE3. As a result we obtain 
the following formula for the exchange energy: 

where So(B) = S(B; [3 = 0). Thus the total ground state energy within our approach reads: 

8(B) =£ (B) + A£(B). (55) 

and the average energy per particle is e(B) = £(B)/N. 

In Fig|| we drew the magnetic field evolution of the average ground state energy per electron e(B) = £(B)/N. The 
three frames correspond to the parameter (3 equal 0.0 (no spin-orbit coupling), 0.3 and 0.6. We find the qualitative 
agreement between our curves and the data reported by Ashoori et ah, obtained in the single electron capacitance 
spectroscopy (SECS) experiments. Comparing the curves in the three frames one can conclude that including the 
spin-orbit interaction brings the model curves fairly close to the measured behavior (the fitting is particularly good 
for/? = 0.3). 



C. Far infra-red absorption 



Let us now discuss the FIR absorption under the magnetic field. Since the magnetic field does not affect the 
structure of the HF wavefunctions, the selection rules ( |io|) remain unchanged and the transition energies are: 



£ f {B) - £\B) = e± = Ml' ± 



(56) 



and we deal with four resonance branches. 

In Fig.^| we compare the dependence of the FIR resonance energies obtained within our model with that reported 
by Demel et alu. Assuming that the experimentally observed higher mode is due to the spin-orbit interaction as 
presented here, we again managed to find a good agreement for (3 — 0.3. Particularly, the zero-field lower resonance 
energy (~ 2.8 meV), the magnetic field at which the anti-crossing occurs (~ 1 Tesla), and the energy at this crossing 
(~ 3.9 meV), seem to be fit very well. A gap separating the anti-crossing levels, observed in the experiment, is 
probably due to a slight anisotropy of the confining potential. 



V. CONCLUSIONS 



The self-consistent theory of a many-electron quantum dot has been developed with the inclusion of the electron- 
electron (Coulomb) and spin-orbit interactions. Included here quantum corrections to the ground state energy of the 
system improve the classical result by Shikin et allEl and compare well with the experiments by Demel et al.u and 
Ashoori et all The spin-orbit interaction is shown to have a strong effect on the electronic structure of the system, 
leading to the splitting of the resonance energy in the far infra-red (FIR) absorption. Predicted here anti-crossing of 
the FIR absorption modes in the magnetic field describes very well the similar behavior reported by Demel et aJj. 

The critical magnetic field, at which the FIR modes cross and a bump in the ground state energy occurs, is around 
1 T, which agrees with the experiments for GaAsCru, and on the other hand with the numerical calculations for a 
two-electron InSb dot by Darnhofer and RosslerEj. 
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APPENDIX A: HARTREE ENERGY AT ZERO MAGNETIC FIELD 



In order to perform the summations over the occupied states in Eqs.([25|-|26|), it is convenient to introduce the non- 
zero temperature of the electrons T, and eventually find the limit for T — > 0. The temperature leads to a replacement 
of the sharp Heaviside function in Eqs.(p5|-p6|) by the smooth Fermi distribution function: 

/ \ ft e + n + +E-TI-— fl s \ 1 /a -,n 

n s {e + n+ + e_n_) = I l + exp + + ^ f=- j , (Al) 

where fi s = // — ^ (e + +6-) and ll is the chemical potential. 

In order to conveniently hide for the moment the constant terms in Eqs.(|2^) and ( |25| ) we introduce the primed 
energy symbol: 

3nNe 2 SnNe 2 9-irNe 2 
E' = E-N- ^— + = E + (A2) 

4eR lOeB, 20eR V ; 

and rewrite Eqs.(p5|-|26|) in the form: 

N = n s (e+n+ + £_n_), (A3) 



nj-Ti — a 



£'{T)= {e+n + + e-n-)-n s {e + n + +e-n-)+N £++£ . (A4) 



£+ +£- 

£+n+ -f- E-ii-) ■ n s ye+n+ -\- E-ji-) -t- is 
n+n—G 

In order to find the thermodynamically stable state we shall further introduce the following thermodynamical potential: 

$= 4>s{e+n + + e-n-), (A5) 



where we define: 4> s (s) = — fcsTln(l + exp Jf ). At low temperatures the minimization procedure with respect to $ 
is equivalent to finding the ground state, as we have: 

£'{T) = $ + (jlN -TS (A6) 

(with the entropy S = d<$>/dT) and: 

£' = $o + e F N, (A7) 

where: $o = h m T— >o $ and ep = hniT— ►() P- 

Let us now introduce the Laplace transformations: 



2 r-C+iOG />0O 

n s{£) = ^— / dpn s (p)e pe , h s (p) = / de n s (£)e" ep , 



Using these expressions we can rewrite Eqs.( A3]- A3) as follows: 



(A9) 



where 



* = £/_ <fe*( e )(-^), (A10) 



2 i-c+ioo £ ep 

2KiJ c _ ioo P(l~e P £ +)(l-e 

1 pc+ioo ep 

zx{e) = — dp— — — ——. (A12) 

2mJ c _ ioo P (1 - e p£ +)(l - e 







We chose the constant c in such a way that all the singular points of the subintegral lie on the left-hand side of the 
integral contour (the contour encloses all singularities inside). Since —dn s /de tends to the Dirac's delta for T — > 0, 
one finds that: 



N = 



i4. 

e+S- 



Pi — + Pi — 



MOs 



£_ 



2 £ - 



Mo. 



«l = l-35 



1 1 

£+ £_ 

_ £+ 
£_ 

2ft 



1 

12 



MPs MPs 

£ + ' £_ 



(A13) 



where /ig s = £f — + an d: 



$ °=4e{^ 



j4 

£ + £- 



+ 



£+ft 



I4)i 



-£_P 2 
MPs Mos 

3 

£ + £_ 



Mas 

£_ 
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-ft 



MOs 



-ft 



Mas 

£_ 



In the above formulae we have introduced the following notation: 

cos(27ria;) 



(A14) 



ft,w = E 



oo 



p 2j -+i(z) = E 



2 2j-lj2 j7r2j ' 

sin(27riir) 
22^27+1^27+1 



Fi(x,y) 



£+£_ 



E 



cos(27rix) — cos(27rjy) 



(«£-) 2 - O+) 2 



00 '|±, sin(27ria;) 



-E 

2J = 1 



2^" 



sin(27rjy) 



(Z£_)2 - (j£ + )2 



(A15) 



The above functions are periodic with the period 1 . In the relevant range of domain the absolute value of each function 
is less then unity. Moreover, note that for < x < 1 



P 1 ( x )~-x+- 



ft(z) 



x 
12 



\Pi(x)\ < l> 
\P 3 (x)\ < 0.009. 



(A16) 



The energy of the system £q one can find from Eq.(Al4) usi ng £f (N) determined from Eq.( A13] ). Taking into account 
that e± = £1' ± aa and introducing: vp = ef / '^/£+£-, Eq.( A13 ) can be rewritten as follows: 



M , 1 (?lO^) 2 +(ao/2) 2 _ 2 

JV "+" 2 (nnj,>-(a /2) 2 _ V F ' 1 



)£+ 



ft 



\ £ + J 12 



(A17) 



The absolute values of all periodic functions are ranged here by unity. Hence, for N ^> 1, we note that all the terms 
under the summation on the right-hand side of the above equation are small compared to N if only: 



KQ,' + a/2 

nn' - a/2 



ft 



MOs 



hn' +a/2J 12 



1 



< 1 or : a < -hQ'r 



(A18) 
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Provided this condition one can use the perturbation method and look for the solution in the form of the following 
series: 



(A19) 



where: \vqf\ \vif\ ^ |^2f| 3> • • • From Eq.( A17 ) we find: 



vqf 



1 (k%) 2 + (a/2) 



2 (hQ' ) 2 - (a/2) 



N 



and, using: /ig s = e 0F — hfl' — v F^(hfl' ) 2 — (a/2) 2 — KQ' : 



V\F 



Pi 



1 

12 



P? 



VPs 
£_ 



Using the above expressions for Vf we find that (for N 3> 1) 



ejr = voF^Jm'of - (a/2) 2 



^1F 



0(iV 



1 

12 



-3/21 



1 + 3^ + 0(^-3/2) 
VQF 



£' = $ + ^£f - y^/«H^[l + 0(i\T 3 / 2 )]. 



(A20) 



(A21) 

(A22) 
(A23) 
(A24) 



The main non-oscillating term in the above formula for $o is of the order of JVy N while the oscillating terms are at 
most of the order of V^V- In the formula for the non-oscillating term is of the order of N^/N and the oscillating 
terms are at most of the order of unity. 

The estimation of the magnitude of spin-orbit coupling constant a is given in Appendix |^. Using Eq. ( ]B5|), the 
final form of the formula for the Hartree energy, given explicitely in Eq. (|27|) , can be now obtained from Eq.(A24) by 
substituting Eq.(B5) and shifting E' back by the previously omitted constant (according to Eq.| 



APPENDIX B: SPIN-ORBIT COUPLING CONSTANT AT ZERO MAGNETIC FIELD 



According to Eq.(g) one needs to estimate the average self-consistent field < SU > acting on an electron. The 

\m*n 12 



energy of a classical particle in the two-dimensional potential ^m*Vl' 2 r 2 , moving on the orbit of radius r, is: 



SU = -m*ttfir 2 + -tt' l. (Bl) 

Substituting the classical variables by the respective operators and averaging SU first over the quantum state [n, m, a}: 

< nma\SU\nma >= -KQ' (n + m + 1) (B2) 
and then over all occupied states, one can obtain the following formula: 

< SU >= ^jWq ^2(n + m+l)- n s (nh%), (B3) 

nma 

where the spin-orbit energy under the distribution function n s has been neglected. One can now use the following 
relation: 

/<9$n\ /<9<f>n\ 

^(n + m + l)-n s (nx + my)= + ( "If- ) ' ( B4 ) 
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and finally arrive at the formula: 



Based on the above equation, the second criterion in Eq.( A18 ) can be expressed as: 



(B5) 



(B6) 



which gives N < 144 for (3 = 0.3. 



APPENDIX C: FOCK ENERGY AT ZERO MAGNETIC FIELD 

Eq.(fj6|) can be written as: 

A £ = -— 'ltT'n 0s (nhn Q )n 0s (n'hn ), (CI) 



e 

n,n' m,rn' 



where the prime in the second sum excludes terms with m = m! , and: 



/" Z = I dr I dr' 



|r — r 

1 f , f , i f dd _j( m _ m ')jflimi( I )Sji l ra'( I )finm( I ')^n'ra'( 3 '') 



, f dxx / dx'x / - -(m-m-Jg ^m^y p-m-v^^mv^y^- m-y^ y^ (C2) 

'o Jo Jo Jo 2n \/ x 1 - 2xx' cos 9 + x' 2 

Expanding the subintcgral into the Legendre polynomials and integrating over 8 one obtains: 



^n,m ;/ / dxX Rum (,x)Rn'm f (^) 



JO 



X 



\ dttR nm (tx)R n/m ,(tx)t lAml jrt 2k a k+lAml a k 
{ Jo k=0 

/do 00 ~) 

dtR nm (tx)R n > m > (tx)t~\ Am \ ^2t~ 2k a k+ \ Am \a k \ 
k=o J 



(C3) 



where Am = m — m! and a k — ^ 2k ) \\ " ■ The integral /" ™ decays rapidly for Am — ► oo and hence we can cut 



off the summation over ml in Eq.(Cl) keeping only the terms with m' = m ± 1. In the integral over x we deal 



with a subintegral with factor e x and thus the most contributes the region around x = 1. Provided that the 



function in braces (in Eq.C3) is smooth with respect to x, it can be replaced by its value at x = 1. Moreover, since 
Rnm(t)Rn'm'(t) ~ e _t the second term in the braces can be neglected as small in comparison with the first one. 
Hence, for |Am| = 1 the integral I^m attains the form: 

II % ~ - / dxx 2 R nm {x)R n , m ,{x) / dtt 2 R nm (t)R n , ml (t) - + —t 2 + . 



L o \Jo 

where we have taken into account the relation 



'o jo Jo V ^ lb 

1 

37 



Iff " \ 2 

dxx 2 R nm {x)R nlm i{x)) , (C4) 



1 2 /'OO 

dtt 2 R nm {t)R n , m ,{t)~- dtt 2 R nm {t)R n , m ,{t). (C5) 
o J Jo 



Since: 



12 



we have: 



and: 



f 

Jo 



dx X 2 Rnm{x)Rn'm±\[x) = ±- r 



n ± m 
1 H 2 — "n',n+l 



^n' ,n- 



n ± m . . 

1 H - I <V,n+l 



n =F to 



A£ — no s (n?ir2o)"-Os("'?i^o) {Inm,n'm-l^m',m-l + ^nm,n'm+l<W,m+l) 

e '0 • • 



nm.n'm' 



4e 2 
3^ 



+ l)no s (n7LQ ) - y - i 



Using Eq. ( B4 ) we finally find the correction due to exchange interaction in the form of Eq. (p 



APPENDIX D: HARTREE ENERGY IN A MAGNETIC FIELD 



Similarly as in the case of the zero magnetic field (see Appendix |A|) it is convenient to shift the Hartree 
£{){B) and introduce: 



E' {B) = £ (B) 



2 D 2 



9irN 2 e 
20ei?(5)' 



given by: 



£q(B) = ( £ + n + + £-n~) ■ n s (e + n + + e_ri_) + N- 



n+n_ <r 



In order to calculate £'^{B) we further introduce the notation: 



v-p = e 



(MT B ) 2 - (a(5)/2) 2 



F ' [{hn' B Y - (a(B)/2) 2 ] 2 - {hw c a{B)/2f 

2 _ [{hn> B f - (a(i?)/2) 2 ] ■ [{hn> B f (^ c /2) 2 ] 

^ p' B ) 2 -( a (B)/2) 2 ] 2 -(ft Wc a(B)/2) 2 



and rewrite Eq.( [Al3| ) as: 



p> 



1 

12 



- 2- 



All the oscillating functions here are small compared to N if: 



Pi 



MPs 

£+ 



1 

12 



P 2 



Mo* 



2-ro , 

e+ e- 



1 

12 



(nn' B ) 2 + {-huj c ) 2 + -hw c 



1 + 2- 



K(B)C7 



<4JV 



1 



1 



(ft^) 2 + (-^ c ) 2 - -fiu, 



1 + 2- 



*(B)t7 
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At low magnetic fields, i.e. for ^huj c <C a(B), the corrections to ft' and a due to the field are negligible and the above 
condition is satisfied for sufficiently high N. At strong magnetic fields, when \Tiuj c ^> a(B), it takes the form: 



2ft', 



2ft', 



< 4A 



2ft', 



2ft', 



which leads to: 



V2ft' B 



< A". 



(D8) 



(D9) 



For example, taking the parameters for GaAs and N = 15 the above yields: -B <C 17 T. The solution of Eq.(D5) can 
be represented as the series: vf = vqf + v\f + • ■ ., where: 



The first order correction we find as: 



-2 



Pi 



11° 



(D10) 



Pi 



Pi 



f'Of 



1 

12 



2F Q 



- 2- 



Pi 



/4 



1 

12 



where [1q s = Eqf — ^ft- Using the above expressions for vp, we obtain: 



sf — Vqf\ 



' [m' B Y - (a(B)/2) 2 ] 2 - (hcu c a(B)/2y 
{hn> B y - {a{B)/2f 



1 + YlL + 0(N~ 3 / 2 ) 
vqf 



£' (B) = -Nv 0F] 



{hn' B y - {a{B)/2 f 



[l + 0(A^- 3/2 )]. 



(Dll) 



(D12) 



(D13) 



Following the similar procedure as for the zero magnetic field (see Appendix |b|) we shall now estimate spin-orbit 
coupling constant, which is now a function of the field. Analogously to Eq.(B3) we have (neglecting here the spin-orbit 
energy): 



< su >= —nn' R • ^ 
2N B n 



(n + m + 1) ■ n s (nMl H — mhuj c ), 



(D14) 



and further: 



a(B) = -i3f B nn' B VN, 



with the renormalizing function: 



(D15) 



(D16) 



z—uj c /2Q' b 



At low fields the function fg tends to unity, while for B — > oo it decays to zero: fg ~ 1/yJV. Therefore, to a good 
approximation, in the definition of fg we can replace ft' B by fto, defined by Eq.(|32|). 
Finally, using the above formula for a(B), and introducing the following function ug\ 



u B = i + {uj c /2n y 



1 



1-4- Z 



l-z\N 'I - z) 
the Hartree energy of the system can be written in the form of Eq. (f48| ) 



(D17) 



z=p 2 P B N/m 
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FIG. 1. The average ground state energy per electron as a function of the number of electrons in the dot. The classical 
result (stars) is taken from Ref. ^l], the experimental data (stars) - from Ref. ^, and the two curves in between (circles) are 
obtained within our model for two values of the spin-orbit coupling constant (3 (GaAs, hu)o = 5.4 meV). 

FIG. 2. The dot radius as a function of the magnetic field for TV = 30 electrons. The three curves correspond to the 
spin-orbit coupling constant /3 equal to 0.0, 0.3 and 0.5 (GaAs, Hujq = 5.4 meV). 

FIG. 3. The average ground state energy per electron as a function of the magnetic field and the number of electrons. The 
three frames correspond to the spin-orbit coupling constant (3 equal to 0.0, 0.3 and 0.5 (GaAs, Tiloq = 5.4 meV). Insets show 
the chemical potentials. 

FIG. 4. FIR absorption spectra of a quantum dot containing 25 electrons. Stars - experiment by Demel et ai.0, lines - the 
model (GaAs, /3 = 0.3, Tiuio — 7.5 meV). 
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